Freezing distributed entanglement in spin chains 
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We show how to freeze distributed entanglement that has been created from the natural dynamics of spin chain 
systems. The technique that we propose simply requires single-qubit operations and isolates the entanglement 
in specific qubits at the ends of branches. Such frozen entanglement provides a useful resource, for example 
for teleportation or distributed quantum processing. The scheme can be applied to a wide range of systems - 
including actual spin systems and alternative qubit embodiments in strings of quantum dots, molecules or atoms. 
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Over the past few years there has been significant interest in 
the propagation of quantum information through spin chains, 
which has been stimulated by potential realizations in sys- 
tems such as large molecules, quantum dot arrays and trapped 
atoms. It is fundamentally interesting to determine just what 
the natural dynamics of these systems permits - and from a 
practical point of view the controlled propagation of quantum 
information is likely to be an essential part of any quantum 
computing device. For macroscopic distances quantum states 
of light are almost certain to form the best medium for quan- 
tum communication - but for shorter, microscopic, distances 
it is equally likely that other media, such as spin chains, could 
make a very useful contribution. For example, in future solid 
state devices there could be a need to provide quantum com- 
munication links over microscopic distances, between sepa- 
rate quantum processors or registers, or between processors 
and memory, analogous to the conventional communication 
that goes on within the computer chips of today. 

Distributed entangled states can be used in order to improve 
the fidelity of a quantum state transmission, which then pro- 
ceeds through teleportation. The entanglement can be pre- 
pared off line, and then purified |1] prior to its application. 
A distributed resource can also facilitate quantum repetition 
12], quantum gates between separated systems or the construc- 
tion of extended entangled resources such as cluster states [3]. 
Previous work has focussed on the production and distribu- 
tion of entanglement but, in order to use it as a resource, it 
is essential that the entanglement can be stored or frozen for 
future operations. In this Letter we present a technique for do- 
ing exactly this, using the natural dynamics of spin chain sys- 
tems. We will use 'spin chain' to cover any physical system 
whose dynamics can be predicted using a Hamiltonian iso- 
morphic to that of a coupled spin chain. This could be strings 
of actual spins (produced chemically, or fabricated) connected 
through interactions, or it could be a string of quantum dots or 
molecules (like fullerenes), containing exciton or spin qubits. 
Strings of trapped atoms form another possibility. 

Production of a useful, frozen, entangled resource requires 
production and good fidelity distribution of entanglement, fol- 
lowed by intervention to freeze it. There has been significant 
study of the propagation of quantum states through spin chains 



or networks. Originally it was shown that a single spin qubit 
could transfer with decent fidelity along a constant nearest- 
neighbour exchange-coupled chain 10]. The fidelity can ap- 
proach unity if the qubit is encoded into a packet of spins 
I.5.] . Perfect state transfer can also be achieved in more com- 
plicated systems, with different geometry or chosen unequal 
couplings 10, 13] ■ Parallel spins chains also enable perfect 
transfer jHI], as does a chain used as a wire, with controlled 
couplings at the ends [9]. Related studies have also been 
made on chains of quantum dots i fioi II], chains of quan- 
tum oscillators ]12] and spin chains connected through long 
range magnetic dipole interactions ]13]. State transfer and 
operations through spin chains using adiabatic dark passage 
has also been proposed l,14J . We will discuss the creation 
of high fidelity spatially separated entanglement — in a form 
suitable for freezing — from the natural dynamics of branched 
spin chain systems, and then show how it can be frozen. The 
consequences of branching were first studied in the context 
of divided bosonic chains, composed of coupled harmonic os- 
cillators [15, 16] - and have since been studied in systems of 
propagating electrons iITtIi . In both of these cases, the dy- 
namics of Gaussian wave packet type excitations have been 
investigated. Our work deals with the propagation of a single 
spin-down excitation localized on a single site, and how this 
moves through a network of spin-up states 

To introduce our formalism, we first consider a one- 
dimensional chain of N spins, each coupled to their nearest 
neighbours. The Hamiltonian for the system is 
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where a\ is the z Pauli spin matrix for the spin at site i and 
similarly (t!|_ = ± icr*. For actual spins, Ei/2 is the local 
magnetic field (in the z-direction) at site i and is the 

local XY coupling strength between neighbouring sites i and 
j + 1. For a coupled chain of quantum dots where each qubit 
is represented by the presence or absence of a ground state 
exciton, Ei is the exciton energy and Ji,i+i is the Forster cou- 
pling between dots i and i + 1 ifioll . The spin chain formal- 
ism applies to both such physical systems, and to any others 
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which can be described by an isomorphic Hamiltonian. The 
computational basis notation for the spin states at each site 
is |0)i = I tz)i and = I iz)i- The total z-component 
of spin (magnetization), or total exciton number, is a constant 
of motion as it commutes with H. It is therefore instructive 
to consider a state of the system consisting of the ground state 
with the addition of a single flipped spin. This state is straight- 
forward to prepare, assuming local control over a spin at, say, 
the end of a chain, and the flipped spin can be regarding as a 
(conserved) travelling qubit as it moves around under the ac- 
tion of the chain dynamics [4j\. An efficient way of represent- 
ing such states in an N spin network is the site basis defined as 
|k) = |0i,02, ...,Ofe_i, lfc,Ofc+i, ...,OAf). A system prepared 
in this subspace remains in it. Now the detailed dynamics de- 
pend on the local magnetic fields or exciton energies, but if 
these are independent of location i then the dynamics favour 
quantum state transfer processes, as akeady mentioned. We 
adopt this limit for our work here. 

The simplest system we consider is a Y structure — used to 
prepare bi-partite entanglement from a simple initial state — 
labelled as shown in the ten site example of Fig. [T] There is 
a deliberate asymmetry in the coupling of the outside sites to 
the hub of the Y. We first focus on the smallest Y-example, 
which has only four sites: one 'input' site, coupled to a 'hub' 
site with strength Ji, and two 'output' sites, coupled to the 
hub with equal strength J2. |— ) = 2^2 (|3) — |4)) is an 
eigenstate of the system, so if we initialize in the state |1) 
then |— ) is decoupled and plays no part in the dynamics. 
We define the orthogonal state |+) = 2"^(|3) + |4)). The 
Hamiltonian in the {|1) , |2) , |+)} space is then: iJ^^ = 
Ji |1) (2| + J2\/2|2) (+|+-ff.c.. The system is thus equivalent 
to a ID three-site spin chain where the 'output' site is the sym- 
metric entangled state 2~2 (|0)3|1)4 + |1)3|0)4). This effects 
perfect state transfer from the input site 1 to the entangled out- 
put so long as the couplings are equal - i.e. if Ji = J2 V2- 

This entanglement creation and distribution extends to spin 
chain systems with longer arms. Using the method of analysis 
of Christandl et al. [Ti], perfect transfer along a longer chain 
is possible through the use of TV — 1 unequal couplings along 
an N site chain, which satisfy 

J,,,+i = a^i{N -i). (2) 

where a is a constant that sets the size of the interaction and 
Ji^i+i is the coupling between site i and site i + 1. As detailed 
in Ref. |01, an array of spins with equal couplings j can also 
be projected onto a ID chain that has couplings Ji , J2 \/2, ex- 
actly as before. According to the general properties of a three 
site linear chain [7], this structure allows for perfect transfer 
between the two extremes so long as Ji = Ji\f2. In this 
case, the perfect transfer of an excitation in node 1 is allowed 
across the Y-shaped structure to the sites at the extremes 3 and 
4, where it must be shared between the two ends, giving rise 
to entanglement. The same \/2 hub branching factor prevents 
hub reflection for propagating wave packets in divided chains 
of harmonic oscillators 1 16], as well as Gaussian wave packets 
of propagating electrons or magnons l,17jl . 




FIG. 1 : (Color online) (a) Ten site Y spin chain network. Its equiva- 
lent one-dimensional representation is shown in (b). 



Extending to larger structures, we define a Y-structure by 
1 1 ^2 , ^3 ) with 1 1 being the number of sites in the input branch 
and the total site number = + Z2 + /a + 1 . As an ex- 
ample, we have performed dynamical simulations using the 
Hamiltonian of the condition and the branching rule 
for the couplings at the hub spin. Fig. [T] shows how the 
hub branching rules applies to the (3, 3, 3) structure, where 
Ji = Jq = aV6, J2 — J5 ^ a^/TO, J3 = a^/T2, J4 = 
for perfect state transfer Fig. |2] shows the result of the sim- 
ulations for both the (3, 3, 3), TV 10, and (10, 10, 10), 
N = 31, structures. The initial condition is ci = 1 for the in- 
put spin (and all others zero), where Cj indicates the amplitude 
coefficient for the state The upper panel of Fig. |2]shows 
the temporal evolution of |cip and the output spin \cn\'^ for 
the two systems. We underline that, due to symmetry, both 
output spins have the same dynamics in this case, and our re- 
sults show that the excitation is completely transferred from 
site 1 to the output sites at time Tr/2a and periodically re- 
turns there at regular time intervals of ir/a. With respect to 
the rescaled time at, the peak corresponding to each revival 
narrows as the length Ik of the branches increases. This is 
important when considering an optimal practical structure for 
distributing entanglement, since a narrower peak puts greater 
time constraints on any entanglement extraction protocol. The 
middle panel of Fig. |2] shows the corresponding fidelity ll20ll 
with respect to the state |+) for the two output spins. Note 
that maximally entangled states can also be produced if the 
length of the output branches is different from the length of 
the input branch, as long as the couplings satisfy Q and the 
hub branching rule. Indeed, in the limit the length of the input 
branch can be reduced to zero, so the excitation is effectively 
made at the centre of an odd N chain. This limit may be help- 
ful for some physical realizations, although having an actual 
input branch may be more practical in some cases. 

We can quantify the entanglement at the output using the 
entanglement of formation, Ep, which measures the number 
of Bell states required to create the state of interest. For a 

two qubit state it is given by: Ep{p) — h ( ^"'"^^'"^ j where 
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FIG. 2: (Color online) Results of simulations on a (3,3,3) and a 
(10,10,10) Y spin chain system, with nearest neighbour couplings 
chosen to satisfy (O and the hub branching rule. Initial condition: 
ci = 1 Upper panel: |cip and |cjvp with respect to the rescaled 
time at. Middle panel: corresponding fidelity with respect to the 
output state |+). Lower panel: entanglement of formation of chain 
ends state. 



h{x) — — xlog2(a;) — (1 — x) log2(l — x) is the Shannon en- 
tropy function, r is the "tangle" or "concurrence" squared: 
T = = [max{Ai — A2 — A3 — A4, 0}]^. The A's are 
the square roots of the eigenvalues, in decreasing order, of 



the matrix pp — p (jy ® '^y P*'^y ® '^y ' where p* de- 
notes the complex co njug ation of p in the computational basis 
|00) , |01) , |10) , 111) U^. Ef between qubits = 7 and 
~ 10 of the (3, 3, 3) structure, and n2 = 21 and ~ 31 
of the (10, 10, 10) structure are shown in the lower panel of 
Fig.|2] as a function of time following initiaUzation of the ex- 
citation on site 1 . As expected, a maximally entangled state is 
obtained after a time 7r/2a and at intervals of 7r/a thereafter. 

To provide a resource, having generated and distributed en- 
tanglement using a spin chain system, it should be isolated or 
frozen in some way. Clearly, as illustrated in Fig.|2] the arrival 
time at the output relative to the preparation time at the input 
is known. At this point, particularly if a number of such states 
are to be purified before use, the dynamics have to be halted. 
This could be done by extracting an entangled state (e.g. with 
swap operations) from the spin chain system and transferring 
it to a storage system of qubits, or to members of spatially sep- 
arated quantum registers or processors, (see e.g. Ol ill ). An- 
other possibility is to physically isolate the two nodes at the 
end of the chain, after the entanglement has formed there. If 
fast local control of the couplings is possible, then the state 
could indeed be frozen at the ends of the output branches. 
To achieve this the couplings have to be switched, fast on a 
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FIG. 3: (Color online) Dynamics of |c,ij (for the output branch 
end spins) versus the rescaled time at, for the bifurcation structure 
shown. This follows the initialization condition ci = 1. 



timescale set by the chain dynamics illustrated in Fig. |2] si- 
multaneously on both output branches. Both of these methods 
present obvious difficulties, and are rather cumbersome, so 
we now discuss a new approach to freezing separated entan- 
glement, which could be very effective for some realizations 
of the spin chain systems. 

A very simple action that can be taken at the end-chain en- 
tanglement arrival time is to apply a phase flip to just one of 
the two output spins, transforming the state to |— ). For ex- 
tended chains this is not an eigenstate. However, from the 
symmetry of the Hamiltonian, its subsequent dynamics can- 
not involve the hub or the input branch, so the subsequent 're- 
vival time' of the state |— ) at the output is reduced. This in 
itself might be useful, but if the structure is modified slightly, 
coordinated action can provide complete freezing. 

The required device has a Y structure, with the ultimate 
spin at each output branch replaced by a bifurcation into two 
final spins. These both couple to the penultimate spin with a 
strength reduced by 1/^/2, compared to the single coupling 
that they replace. In terms of the four end spins of the full 
device, the natural dynamics will generate a state of the form 
IV,) = i(|0, 0,0,1) + 10,0,1,0) + 10,1,0,0) + |1,0,0,0)). 
The structure and dynamics are shown in Fig. [3j where we 
plot the branch end spin probabilities |c„. p as a function 
of rescaled time at. If at a probability peak a phase flip 
is applied to one spin out of each pair, a state of the form 
iVa) = i(|0, 0, 0,1) - |0, 0,1, 0) + |0,l, 0,0) -|1, 0,0,0)) 
results. This is an eigenstate of the system and the entangle- 
ment is thus frozen. Although four spins are involved, the 
spatial separation is between two pairs of spins. Each pair can 
be viewed as a storage buffer for one qubit. The system con- 
tains spatially separated and stored bipartite entanglement, 
which could be released for future use by single qubit opera- 
tions and/or coupling to other systems at the branch ends. 

To further quantify the entanglement freezing, we investi- 
gate the dependence on the timing of the coordinated phase 
flips. The upper panel of Fig. |4] shows the fidelity of the state 
as a function of time, measured against the desired frozen state 
\4>a), for five different choices of the (simultaneous) phase flip 
timing. This shows that even if the timing is not perfect a high 
proportion of the state is frozen. The middle panel shows the 
subsequent dynamics of the remains of the state, for the same 
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FIG. 4: (Color online) Behaviour of the frozen entanglement as a 
function of phase flip timing, for the system of Fig. [3] (a) Time 
dependence of the fidelity against \^a} for five different flip times, 
assumed to be simultaneous with each other, (b) Time dependence 
of the fidelity against \i}}s) for the same five flip times, (c) Frozen 
bipartite entanglement of formation Ef as a function of the times ti 
and t2 at which the phase flips are applied. 



five examples of phase flip timing, through plots of the fidelity 
against tiie un-flipped state lips)- These amplitudes continue 
to propagate periodically through the spin chain system, with 
an overall normalization set by the timing of the phase flips. 
The lower panel quantifies the amount of bipartite entangle- 
ment frozen, showing the entanglement of formation Ep as a 
function of the times ti and t2 of each pulse. The encoded ba- 
sis used for these calculations uses states of the pairs of spins 
at each chain end and takes the following form: |0l) = |00) 
and|lL> = 2-i/2(|01) - |10)). Clearly if the phase flips are 
fast compared to the system dynamics, and can be timed suffi- 
ciently accurately and simultaneously compared to the period 
of these dynamics, it is possible to freeze a high quality spa- 
tially separated Bell state, encoded in the two pairs of spins. 

To conclude, we have shown how entanglement — created 
and distributed in branching spin chain systems — can be 
frozen, encoded into pairs of spins located at the ends of Y 
networks. Such distributed entanglement provides a useful re- 



source, for example for teleportation or distributed quantum 
processing. In contrast to the use of spin chains to propagate 
quantum states from one place to another with as high a fi- 
delity as possible, we see some advantage in building up a 
high fidelity entangled resource off line. Real systems, with 
their inevitable imperfections, will almost certainly degrade 
transmission fidelities, even if in principle these approach 
unity. Certainly, with the 'off-line' resource approach, pu- 
rification in could be applied to build up a higher fidelity 
resource than can be achieved by direct transmission. This 
can then be used to transfer quantum states or some form of 
quantum communication. In effect, the concept of a quan- 
tum repeater |2] could be employed in a solid state, spin chain 
scenario. The key ingredient — entanglement distribution and 
freezing — ^results from the basic dynamics of the branched 
spin chain systems, simply prepared initial states and coordi- 
nated single qubit operations. This opens new possibilities for 
quantum communication in solid state systems, especially as 
fabrication or creation of solid state systems that can operate 
as spin chains continues to progress. 

B. W. L. acknowledges support from the Royal Society and 
QIPIRC (No. GR/S82176/01). 



* Electronic address: Iida500@york.ac.uk| 
[1] C. H. Bennett, et al., Phys. Rev. Lett. 76, 722 (1996). 
[2] H.-J. Briegel, et al., Phys. Rev. Lett. 81, 5932 (1998). 
[3] R. Raussendorf and H. J. Briegel, Phys. Rev. Lett. 86, 5188 
(2001). 

[4] S. Bose, Phys. Rev. Lett. 91, 207901 (2003). 

[5] T. J. Osborne and N. Linden, Phys. Rev A 69, 052315 (2004). 

[6] M. Christandl, et al., Phys. Rev. Lett. 92, 187902 (2004); G. 

De Chiara, et al., Phys. Rev A 72 012323 (2005).M.-H. Yung 

and S. Bose, Phys. Rev A 71 032310 (2005); R Karbach and J. 

Stolze, Phys. Rev. A 72, 030301(R) (2005). 
[7] M. Christandl, et al., Phys. Rev A 71, 032312 (2005). 
[8] D. Burgarth and S. Bose, Phys. Rev A 71, 052315 (2005); D. 

Burgarth and S. Bose, New J. Phys 7, 135 (2005). 
[9] A. Wojcik, et al., Phys. Rev A 72, 034303 (2005). 
[10] I. D'Amico in Focus on Semiconductor Research (edited by 

Nova Science Publishers, 2006), (Preprint cond-mat/051 1470 1. 
[11] T. P Spiller, I. D'Amico and B. W. Lovett, New J. Phys. 9, 20 

(2007). 

[12] M. B. Plenio and FL. Semiao, New J. Phys. 6, 36 (2005); M. J. 
Hartmann, M. E. Renter and M. B. Plenio, New J. Phys. 8, 94 

(2005) . 

[13] M. Avellino, A. J. Fisher and S. Bose, Phys. Rev A 74, 012321 

(2006) . 

[14] A. D. Greentree, S. J. Devitt and L. C. L. HoUenberg, Phys. 

Rev A 73, 032319 (2006); T. Ohshima, et al.. Preprint 

quant-ph/0702019 
[15] M. B. Plenio, J. Harfley and J. Eisert, New J. Phys. 7, 73 (2004). 
[16] A. Perales and M. B. Plenio, / Opt. B 7, S601 (2005). 
[17] S. Yang, Z. Song and C. P Sun, Eur. Phys. J. B 52, 377 (2006); 

S. Yang, Z. Song and C. P Sun, Front. Phys. China 1, 1 (2007). 
[18] S. J. Devitt, A. D. Greentree and L. C. L. Hollenberg, Preprint 

quant-ph/05 11084 
[19] M. Patemostro, H. McAneney and M. S. Kim, Phys. Rev. Lett. 



5 

94, 070501 (2005). Information p. 409 (Cambridge University Press 2000). 

[20] For our calculations we adopt the fidelity definition of M. A. [21] W. J. Munro, et al, Phys. Rev. A 64, 030302 (2001). 
Nielsen and I. L. Chuang, Quantum Computation and Quantum 



